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ABSTRACT 

The recent discovery of the first four afterglows of short-hard gamma-ray bursts (SHBs), and the 
properties inferred for these bursts - as well as for four other SHBs with known or constrained redshift 
- suggests that these events typically result from long-lived progenitor systems. The most popular 
model invokes the merger of two compact objects, either double neutron star (DNS) binaries or neutron 
star-black hole (NS-BH) systems. Such events should emit a significant fraction of their ultimate 
binding energy in gravitational radiation in the frequency range that is accessible to current and next- 
generation ground-based gravitational-wave observatories. In this work we combine the census of SHB 
observations with refined theoretical analysis to perform a critical evaluation of the compact binary 
model. We then explore the implications for gravitational wave detection of these events. Beginning 
from the measured star-formation rate through cosmic time, we consider what intrinsic luminosity 
and lifetime distributions can reproduce both the known SHB rcdshifts and luminosities as well as 
the peak flux distribution of the large BATSE SHB sample. We find that: (1) The typical progenitor 
lifetime, r*, is long. Assuming lognormal lifetime distribution t, > 4[1] Gyr (at 2[3]cr c.l.). If the 
lifetime distribution is a power-law with index 77 then 77 > — 0.5[— 1] (at 2[3]cr c.l.). This result is 
difficult to reconcile with the properties of the observed galactic DNS population, suggesting that if 
SHBs do result from DNS mergers then the observed galactic binaries do not represent the cosmic 
DNS population. (2) We find that the local rate of SHBs is at least TZshb ^ 10 Gpc~^ yr~^ and 
may be higher by several orders of magnitude, significantly above previous estimates. (3) We find 
that, assuming that SHBs do result from compact binaries, our predictions for the LIGO and VIRGO 
event rates are encouraging: The chance for detection by current facilities is not negligible, while a 
coincident observation of electromagnetic and gravitational radiation from an SHB is guaranteed for 
next-generation observatories. 

Subject headings: gamma rays: bursts; gravitational waves 



1. INTRODUCTION 

More than a decade ago it was realized that Gamma- 
Ray Bursts (GRBs) can be divided into two well-defined 
sub-populations: the majority of the observed bursts 
(about 3/4) had durations longer than ~ 2s and rel- 
atively softer observed spectra, and the minority (~ 
1/4) had short durations and harder observed spectra 
ijKouveliotou et al~\\l99^ . Accordingly, the two compo- 
nents of this bimodal population are often referred to as 
long-soft GRBs and short-hard GRBs (SHBs hereafter). 

Observational evidence accumulated in the last 
few years has conclusively shown that long-soft 
GRBs are associated with supernovae (SNe), 
and result from the death of short-lived massive 
stars ( G alama et al. 1998; Kulkarni et al. 1998; 
Bloom et al1ll999: Bloom. Kulkarni & Diorgovski 2002; 



Stanek et al. Il2nn,3; .Hiorth et al 



iMatiiggon_ 



2nn3l ILmkin .^a/. l l2nn4t IGal-Yarn g/. I M 
Malesani et al. ll2004HCobb et al. ll2004HThomsen et al. I 
2004; Sodcrbe rg et al. I \200g^ . Comparable studies of 
SHBs were not conducted due to observational dif- 
ficulties in obtaining accurate localizations for these 
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events. 

A breakthrough in the study of SHBs occurred ear- 
lier this year when the Swifi^ and HETE-^ satel- 
lites provided the first timely and accurate localiza- 
tions of SHBs 050509b, 050709 and 050724, leading to 
the detection of afterglow emission in X-ray, optical, 
IR and radio wavelengths, as well as to the identifi- 



(IBloom et al l200,5HKulkami et al l2nn,'llGehrels et al 1 




Castro- Tirado et al. 1 




Prochaska et al. ||2005l 



Berger et a l. 2005; Hiorth et al. 2005). These studies 



provide three strong indications that SHBs result from a 
different type of physical progenitor system. First, in all 
three events there are no indications of an associated SN, 
to within strict limits. SNe similar to those associated 
with long-soft GRBs would have been easily detected. 
Second, two of these events occurred within early-type 
galaxies, with little or no recent star formation. None of 
the many tens of long-soft GRBs with comparable data 
are detected in such galaxies. Finally, all of these bursts 
occurred relatively nearby, at redshifts z < 0.3, in con- 

^ http:/ /swift. gsfc.nasa.gov/docs/swift/swiftsc. html 
3 http: //space. mit . edu /HETE /Welcome. html 
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trast with long GRBs that are typically at z > 1. Taken 
together these observational facts strongly indicate that 
SHBs have a different type of progenitor system, and that 

some of these progenitors are long-l ived (> 1 Gyr) . 

The most popular model (e.g. , lEichler et al. 1 1198^ 
iNaravan. Paczvnski fc PirarJll992) for SHBs invokes the 
merger of two compact objects, such as double neutron 
stars (DNS) or a neutron star and a black hole (NS- 
BH). Such models are appealing because they predict 
events with comparable timescale and energy release to 
those observed in short bursts. Such considerations, com- 
bined with the recent results showing that SHBs require 
a long-lived progenitor system, led several groups to sug- 
gest that SHBs result from coinpact binary mergers (e.^ 



iBloom et a l. '2005' iFox et al . "2005*: ' Berger et ani20o! 
More detailed comparisons with model predictions were 
limited by the small numbers of SHBs with known host 
galaxies a nd redsh i fts. 

Gal- Ya m et al. I l)2005D have expanded the sample of 
SHBs available for study using new and archival ob- 
servations of "historical" SHBs, allowing them to es- 
tablish a statistically significant association between 
two additional SHBs (790613 and 000607) and their 
probable host galaxies or clusters. Combining these 
data with observations of the recent Swift and HETE- 
2 SHBs discussed above, and an additional Swift 
burst (SHB 050 813) possibly lo c ated i n a z = 0.72 
galaxv cluster llGladders I 120051 iBergeil 120051 
iProchaska et al. I 120051) . iGal-Yam et al (2005) con- 
fronted predictions from binary NS merger models ijAndol 
11)04; Guctta & Piran 2005) with the properties of the 
observed SHB sample. They found that the model pre- 
dictions are difficult to reconcile with the properties of 
the observed SHB sample, since the range of models con- 
sidered in the literature inevitably predicts a large frac- 
tion of short-lived merging systems, which will preferably 
reside in late- type hosts and occur at higher redshifts, in 
contrast with the observed SHB sample which is domi- 
nated by events located in early-type hosts and at lower 
redshifts. 

In this work we apply a refined theoretical analysis to 
the Gal- Yam et al. I 1 )200 5') sample in order to constrain 
the local rate and the lifetime of SHB progenitors. Our 
goal is to confront our findings with compact binary mod- 
els for SHBs and discuss the implication for gravitational 
wave detection from SHBs by the LIGO and Virgo ob- 
servatories. 

The structure of the paper is as follows: We develop 
the theory that we need for the analysis of the data in 
§2. The observed sample, its analysis and the derived 
constraints are presented in §3. We compare our results 
to the predictions of compact binary mergers in §4 and 
discuss the prospects for detecting GWs from SHBs un- 
der the assumption that they originate from a DNS or 
BH-NS mergers in §5. We summarize our conclusions in 
§6. 

2. THEORY 

We use the combination of the two-dimensional ob- 
served redshift and luminosity distribution together with 
the observed SHB flux distribution (logN-logS), to con- 
strain the intrinsic SHB rate at redshift z, Rshb{z) 
and the intrinsic luminosity function ML). This method 
is an extension of previous works l|Piranl Il992t lAndol 



120041 iGuetta fc PiranI l2005j) that used only the one- 
dimensional observed redshift distribution and the peak 
flux distribution. Although we focus here on SHBs, the 
method is applicable to long-soft GRBs as well. 

The two-dimensional observed redshift and luminos- 
ity distribution is derived from the intrinsic distributions 
via: 



dNo. 



bs 



dLdz 



1 + z dz 



(1) 



where Nobs is the observed SHB rate and </)(-£) is the 
intrinsic peak luminosity function (which we assume to 
be independent of z). < S{P) < 1 is the probability for 
detection (including redshift determination) of a burst 
with a peak photon flux P which in turn depends on L 
and z as well as on the spectrum of the bursts. Rshb{z) 
is the intrinsic SHB rate per unit comoving volume and 
comoving time. Since SHB progenitors are most likely of 
a stellar origin we expect: 



/■°° dt 
Rshb{z)o, j SFR{z')f{t{z)-t{z'))—dz', 



(2) 



where SFR{z) is the star formation rate at redshift z 
(per unit comoving volume and comoving time), t{z) is 
the age of the universe at redshift z, and /(r) is the 
fraction of SHB progenitors that are born with a lifetime 
r. 

S{P) can describe a single detector or a combination 
of several detectors, each weighted by its field of view 
and operational time. In principle if S{P) is well known 
and if the observed sample is large enough then the in- 
trinsic distributions can be extracted from Eq. ^ In 
reality we have to work with a limited sample as well 
as poorly understood S{P). In the case of GRBs (long 
and short) there is a large sample of bursts, observed 
by BATSE'', for which only the peak flux distribution is 
available while the redshift (and thus luminosity) is un- 
known. The BATSE sample can constrain the intrinsic 
distributions by considering the observed flux distribu- 
tion which is an integration of Eq. ^ 



dz 



dL^^ (3) 



where 



dP dP Jq J L,min{z,P) dLdz 



i™n(z,P) -4^difc(z)P. 



(4) 



dL{z) is the luminosity distance and k{z) is It depends on 
the spectrum of the bursts and includes the k-correction 
as well as the conversion from energy flux to photon flux. 
k{z) is assumed to be a function of the redshift only. 
Note that the LogN-LogS distribution is iV(> P) in our 
notation. 

If (j){L) is a single power-law, (j){L) = (j)QL~^ , with no 
upper or lower cutoff (within a luminosity range that we 
discuss below) then the integral over z in Eq. |21does not 
depend on P and thus the observed peak flux distribution 
does not depend on Rshb and simply satisfies: 



dNobsjP) 
dP 



oc P-^S{P). 



(5) 
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Similarly, the integral over L in eq. ^ results in: 



bs 



dz 



= (4^4fc(z))i-^0o 



RSHB dV 
1 + z dz 



J p-^S{P)dP, 



(6) 

thereby eliminating the dependence on the detector 
thresholds. Naturally for /? < 2 an upper cutoff must 
exist while for /? > 1 a lower limit is necessary. However, 
if the lower cutoff is low enough so that it affects only a 
negligible volume, and if the upper cutoff is high enough 
so it affects only the detection at high redshift, then Eqs. 
© and © are applicable (these cut-offs also prevent the 
integral over P in Eq. Elfrom diverging). Therefore, if 
the observed peak flux distribution can be fitted by Eq. 
El then the luminosity function can be a single power- 
law. In this case we can readily use data sets for which 
S{P) is not well known. In such cases Eq. enables 
a comparison of the one-dimensional observed redshift 
distribution with model predictions. Unfortunately the 
observed luminosity distribution depends on S{P) even 
when the luminosity function is a single power-law. If 
S{P) is well known, a better constraint on the intrinsic 
distributions can be obtained by a comparison with the 
two-dimensional luminosity- redshift distribution (Eq. ^ . 

The above formalism is applicable to any astrophysical 
transient as long as its detectability depends only on its 
peak flux. Eqs. n]and|2|are the most general, assuming 
only that the luminosity function does not evolve with 
the redshift. Eqs. and are applicable only when 
the luminosity function is a single power-law and k{z) 
depends only on the redshift. 

3. CONSTRAINTS ON PROGENITOR LIFETIME AND THE 
RATE OF SHBS 

In this section we use the methods presented in fJ21to 
analyze the observed sample of S HBs. We base our stucl y 
on the SHB sample compiled bv lGal-Yam et aEI l)2005D . 
This sample includes 4 relatively recent SHBs localized 
by Swift and HETE-2, and 4 "historical" bursts local- 
ized by the IPN over the years, and studied in retrospect 
by these authors. We summarize the known redshift 
an d host galaxy p r opert ies of each burst, as compiled 
bv lOal-Ya'm 7tlX\ IpMh . in Table 1. 

When using this sample for statistical studies of the ob- 
served redshifts, as we do here, one has to account for se- 
lection biases. Specifically, a selection bias that disfavors 
the identification of the host galaxy of an SHB at high 
redshift may skew our results. Fortunately, the sample 
we use is almost complete, making correction for selec- 
tion biases simple. At the time that the analysis was first 
made Swift had detected and made prompt follow- up ob- 
servations of only four SHBs. Three are included in our 
sample while the last (SHB 050906) mi ght be associated 
with a galaxy at a distance of 140 Mdc l|Levan & Tanvl^ 
|2005() . If correct, this low redshift strengthens our re- 
sults, as we shall see below; however, we do not include 
this burst due to the association being unconfirmed at 
this time ^ . 

^ During the refereeing process Swift localized adittinal handful 
of SHBs. For most of these bursts there is currently no redshift 
or host information. For SHB 051221 the redshift and host are 
known (Sodcrbcrg ct al. 2006). We stress that this burst, or future 
bursts, cannot be added to a sample that is used in a statistical 
comparison to the observations without proper consideration of the 



TABLE 1 

Host galaxies and redshifts of SHBs 



SHB 


Redshift 


Host Galaxy Association 


Ref. 




z 


Type 


significance 




790613 


0.09 


E/SO 


~ 3cr 


1 


000607 


0.14 


Sb 


~ 2(7 


1 


050509b 


0.22 


E/SO 


3-4:(T 


2 


050709 


0.16 


Sb/c 


Secure 


3 


050724 


0.26 


E/SO 


Secure 


4 


050813 


0.72t 


E/SO 




5 


001204 


> 0.25|0.06J 






1 


000607 


> 0.25[0.06] 




l[2]a 


1 


Note. 


— 1. ICial-Yam et 


al. 1 1200,1) 2. IBloom 


et al. 1 



1120051} ; (Kulkarni et al. ' fSOOF); "Castro-Tirado et al. f2003); 
"Oehrels et a/. 12005) 3. Fox et al. 12005) 4. Bcrgcr et oTl 
(2005); Prochaska et al. I 120051') 5. iGladders et al. I (200^ 

Berger (2005); Proch aska et al 1 120051) 

t while the paper was in the referring process, IBergeil 
1 2006)suggested that this burst might be actually associated 
with a cluster at a redshift ft: 1.8 



The 4 IPN SHBs studied by IGal-Yam et dl lEM) 

comprise a complete sub-sample of IPN SHBs local- 
ized to within 10 arcmin^, defined by these authors 
based on a-priori technical issues (e.g., galactic latitude) 
which should not bias the physical properties of t hese 
bursts. Of these 4 bursts. Gal- Ya m et al. I l)2005(l de- 
termine probable redshifts and host galaxy types for two 
events. For the other two events they set a lower redshift 
limit. This limit is based on a null detection of luminosity 
overdensities in the fields of these SHBs, and thus n o host 
galax y information is available. Since [Gal- Yam et al. I 
120051 used a luminosity-based test for associating hosts 
with SHBs, their search may be biased against bursts in 
low-lumiosity galaxies, and against bursts at high red- 
shift. To avoid the bias mentioned above our statisti- 
cal analysis is based on the density distributions derived 
from the lower limits on the redshift of these two bursts 
(sec 3!rT|l . 

We use the observed peak flux sample in the current 
BATSE catalog^. Therefore it is convenient to define 
the flux and intrinsic luminosity quantities discussed in 
f|21 in an energy range that corresponds to the BASTE 
window for detection, 50 — 300 keV. In this window the 
spectrum can be well approximated as a single power 
law F,j oc w here values range between — 1 < a < 
0.5 ijGhirlanda. Ghiscllini & Celotti 2004) with a typical 
value a « —0.5. Using this spectrum we obtain fc(z) « 



2 X 10-'^(1 -|-z)-(i-") 



erg. 



3.1. 



The progenitor lifetime of SHBs 
3.1.1. Single power-law luminosity functions 
Following the discussion in §|21we first try to fit the ob- 
served peak flux distribution "^^y^-* using a single power- 
law with an index /3. From the current BATSE catalog 
we extract all the short (Tgo < 2s) bursts with peak flux 
in the 64 ms timing window of Pq^ > 1.5 ph/cm^/scc (at 
this range S{P) « 1 for BATSE), resulting in a hst of 340 
bursts. We calculate the maximum likelihood and find 

selection effects that prevented the redshift determination of the 
rest of the localized SHBs 

^ http:/ /www. batse.msfc.nasa.gov/batse/grb/catalog/currcnt/ 
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that [3 = 2±0.1 provides the best fit. The x^/d.o.f of this 
modei is 1.17, confirming it is indeed a good fit. Using 
Eq. |S1 we therefore initially consider a single power-law 
luminosity function 4>{L) oc . 

The 8 SHBs with known (or constrained) redshifts in- 
cluded in our sample were detected by several different 
instruments (see Table 1). The detection threshold for 
each experiment is not well known and construction of 
S{P) for this sample is therefore currently impossible. 
Hence, the full sample can be used only for comparison 
with models of the redshift distribution based on a single 
power-law luminosity function. Should a larger sample 
of bursts for which the detector sensitivity is well un- 
derstood become available in the future, the more con- 
straining comparison with model predictions using the 
two-dimensional L-z distribution (Eq. ^ could be car- 
ried out. 

Following ISch^^ lEffll and lCxuetta fc PiranI lEffll 

we base our redshift distribution model calculations on 
the star formation history (SFH) parametrization SF2 of 
IPorciani fc Mad au (2001): 



10 



exp(3.4z) [n„^{l + z)^ + Ha]^^^ 
exp(3.4z) + 22 (1 -|- z)3/2 



(7) 



and we adopt the standard cosmology ( flm = 0.3, 
flA = 0.7 and Hq ^ 70 kms^^ Mpc^^). The probability 
distributions we consider for the lifetime are lognormal. 



/(r)dr = {TaV2^)-^exp[-{ln{T)-ln{n))'^/2a'^]dT with 
various values of and narrow (cr = 0.3) or wide (cr = 1) 
dispersions, as well as power-law distributions /(t) cx r'' 
with a lower cutoff at 20 Myr and an upper cutoff that 
is larger than the Hubble time (our results do not de- 
pend on the exact values of the cutoffs). For each of 
the se models we carry out a maximum likelihood analy- 
sis l|Press et al. 111992) in order to constrain the typical 
delay (r*) or power-law index {rf). We do so by calculat- 
ing the likelihood of the observations (table 1) for every 
set of model predictions (calculated using Eqs. I21E1&IZI)- 
The two bursts with redshift lower limits are accounted 
for in the following way. For each burst we calculate the 
functional form of the lower limit, l(z) which is defined 
as 1 minus the significance level in which the hypothesis 
that the burst is at redshift z is rejected. So, l(z) goes to 
1 for high redshifts (the possibility that these bursts are 
at high-z is not constrained) and l(z) decreases with z, 
for example, assuming the value 0.05 at z=0.06 (table 1). 
The met hod use d to calculate l(z) is described in detail 
2005) and is applied to these bursts in 
5005f'l . The likelihood of a burst with a 
lower limit is the probability that an observed burst will 
be consistent with the observed lower limit for a given 
model of redshift distribution: J {dNobs/dz)l{z)dz. 

The results of the likelihood analysis that we carry 
for the different lifetime distributions are presented as 
one-sided probability curves in Figure ^ For a model 
likelihood £ the one-sided probability is defined as (1 — 
P^2^i{2ln[Cmax /£-]))/ 2 where Cmax is the maximal like- 
lihood and Px^.i is the cumul ative distributio n with 
one degree of ifreedom (e.g., iPress ef "ZIIT99I . Fig- 
ure n shows clearly that long life-times are favored. For 
the narrow lognormal distribution we find that the most 
likely delay is r, — 6.5 Gyr and its 5% — 95% confidence 
interval is 4 < < 9.5 Gyr, while the probability that 
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Fig. 1. — Left: The one-sided probability curves, which result 
from the maximum-likelihood analysis, for a narrow (thick line) 
and wide {thin line) lognormal lifetime distributions. For lifetime 
values, T, that are smaller [larger] than the most probable value, 
the one-sided probability is the probability that r* < r [r* > r]. 
For the most probable value this probability is 0.5. Considering a 
narrow lognormal distribution the most likely delay is t» =6.5 Gyr 
and its 5% — 95% confidence interval is 4 < r, < 9.5 Gyr, while the 
probability that t» > 1 Gyr is 99.98%. Assuming a wide lognormal 
distribution the most likely delay is r« = 20 Gyr while r« > 4[1] 
Gyr at 95 [99. 9]%. Right; The one-sided probability curves of a 
single power-law lifetime distribution. The most probable power- 
law index is jj = 0.6, the 5% — 95% confidence interval i s —0.5 < 
r] < 2.6 and the probability that r; > -1 is 99.5%. See 1.11 for 
details 

r* > 1 Gyr is 99.98%. Assuming a wide lognormal dis- 
tribution we find that the most likely delay is larger than 
the Hubble time (r, = 20 Gyr), that r > 4 Gyr at 95% 
confidence and that the probability that > 1 Gyr is 
99.9%. Considering a power-law lifetime distribution we 
find that the most probable power-law index is 77 = 0.6, 
the 5% — 95% confidence interval is —0.5 < rj < 2.6 and 
the probability that 77 > -1 is 99.5% 

The normalized (mtrmszc) redshift distributions, 
Rshb{z), resulting from various lifetime distributions 
are presented in Figure |21 Comparison between the 
predicted and observed cumulative redshift distributions 
for several representative models are depicted in Figure 
13 This figure vividly illustrate the results from our 
maximum likelihood analysis, namely that models con- 
sistent with the data must have a long typical delay(> 4 
Gyr). Models with no typical lifetime (power-law 
distributions) must have a birthrate of progenitors per 
unit logarithmic lifetime that increases significantly as 
a function of lifetime (i.e. rj > —0.5). This figure also 
shows that indeed all the models we considered provide 
a good fit to the data when t*(77) takes its most likely 
value. It also demonstrates that the reason that models 
dominated by short-lived systems are rejected is that 
they under-predict the fraction of bursts at low redshift 
(z<0.3). 

Given that the sample that we use is small, and to test 

^ In these calculations we took the redshift of SHB 050813 as 0.72 
IGladders et al. l2005H Berger''2005': "Prochaska et al. 2005). while 
the paper was in the referceing process, Bcrecr (2006) suggested 
that this burst might actually be associated with a cluster at a 
redshift higher than 0.72. Repeating the analysis while taking a 
lower limit of z > 0.72 for SHB 050813 does not change the results 
significantly. For example, for narrow lognormal distribution, the 
most likely delay is r, = 5.8 Gyr, the 5% — 95% confidence interval 
is 3.5 < T* < 9 Gyr while the probability that r* > 1 Gyr is 99.8%. 
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Fig. 2. — The normalized intrinsic SHB rate (Rshb) a func- 
tion of the redshift for several lifetime distributions and SFH pa- 
rameterizations (see legend). 

for biases in our maximum likelihood analysis, we have 
carried out a jackknife analysis (Efron 1982). Namely, 
we repeat the analysis 8 times, in each occasion remov- 
ing a different burst from the sample (i.e. using a dif- 
ferent sample of 7 bursts). We find a negligible bias, for 
example the bias of the most likely r, [77] for a narrow 
lognormal [power-law] distribution is —0.05 Gyr [0.15], 
implying that the most likely values that are obtained by 
the maximum likelihood analysis are not biased because 
of the small sample size. We also find that our results are 
not driven by any single burst. For example, for a nar- 
row lognormal lifetime distribution the most likely value 
of T* in the different samples varies between 5.8 Gyr and 
8 Gyr (in the full sample it is 6.5 Gyr) and in all the 
samples t* > 3.4 Gyr at more than 95% confidence level 
while the probability that > 1 Gyr is, in all samples, 
higher than > 99.8%. 

Given the observational uncertainty in the SFH at high 
redshift we repe at our analysis using the SFH formula 
SFl from iPorciani k Madaul 1)200 1)) : 



SFRi{z) OC 



exp(3.4z) [nmii + zf + n. 



11/2 



exp(3.8z) + 45 



(1-f z)3/2 



(8) 



In this model the star formation rate falls exponentially 
at redshifts larger than 1.5. The results we obtain using 
this SFH model are similar to those obtained using Eq. [3 
For example considering a narrow lognormal distribution 
with SFRi we find that the most likely delay is r* = 6 
Gyr and its 5% — 95% confidence interval is 3.5 < < 9 
Gyr while the probability that > 1 Gyr is 99.95%. 
The reason is that the models differ only at high redshift, 
with little impact on the distribution at low redshifts. Of 
the three models presented in IPorciani fc iVladaul i 



this model places the most conservative constraints on 
/(r) as it predicts the largest fraction of SHBs at low 
redshift (for a given lifetime distribution). Therefore, 
lifetime distributions that are ruled out when applied 
to this model will be ruled out with great er significance 
when the other IPorciani fc Madaul l|2001f) SFH models 
are used. 

We return now to the upper and lower cutoffs of the 
luminosity function. As described above, for a given se- 
lection of Rshb and S{P) there are critical values of the 
lower (upper) luminosity function cutoffs, below (above) 
which the observations are insensitive to the value of 




Fig. 3. — The cumulative observed redshift distribution as pre- 
dicted by various lifetime distributions when the luminosity func- 
tion is <f>{L) OC and the star formation history is given by Eq. 
171 For each functional form of the lifetime distribution we present 
the most likely distribution and the distribution that is rejected at 
95% confidence level. The cumulative redshift distribution of the 
observed data (shaded area) is bracketed between the lower solid 
line, which is the cumulative redshift distribution of the six bursts 
with known redshifts, and the upper solid line, which includes also 
the contribution of the two bursts with upper limits, as given by 
l(z) (see text). The figure demonstrates that the most likely models 
provide a good fit to the data while the models that are rejected 
undcr-predict the fraction of bursts at low redshift 2 < 0.3 



the cutoff. Among all the models of Rshb that we 
have examined, and for threshold values that vary be- 
tween* 1 — 50 ph/cm^/s, the fit to the LogN-logS dis- 
tribution, as well as the fit to the observed redshift dis- 
tribution, are insensitive to a lower cutoff Lmiti < 10*^ 
erg/s and to an upper cutoff Lmax > 10^^ — lO^'^ erg/s. 
Namely, a luminosity function with 4>{L) oc for 
10"'^ < L < 10^^ erg/s and 4>{L) = otherwise is con- 
sistent with all available data. The luminosity of every 
SHB in our sample falls within this range as well. Al- 
though the sample that we consider is consistent with 
Lmin — lO'*^ erg/s, recent results suggest a lower Lmin- 
First, if SHB 050906 is associated with t he nearby galaxy 
IC 32 8 at 140 Mpc, as suggested by iLevan fc TanviJ 
|200 5l), then its luminosity is about 10''^*erg/s. Second, 
iiTanvir et al. I ((2005) found evidence that at least 5% of 
the SHBs detected by BATSE are at z < 0.025 (« 100 
Mpc). Interestingly, a comparable fraction (« 5%) is 
predicted by our best fit model {(j){L) cx L"^ and a typ- 
ical lifetime of 6 Gyr) if the lower cutoff is as low as 
Lmin < 10^^ erg/s. A slightly smaller, but consistent, 
fraction (« 3%) is predicted if Lmin = 10^'' erg/s while 
the fraction drops to less than 0.5% for Lmin — 10'** 
erg/s. 

If Lmin ^5 10^^ erg/s then the SHB luminosity func- 
tion overlaps with the energies observed in giant flares 
from soft gamma-ray rep eaters (SGRs) (jHurlev et al] 
120051 iPalmer et 01112005(1 . This overlap in luminosi- 
ties raises (again) the question whether SGRs can pro- 
duce flares with isotropic equivalent observed luminosi- 
ties much greater than ^ 10^^ erg/s, as we previously 
suggested in lNakar et al. 1 1|2005|) . In this case giant flares 
from SGRs can be the source of all observed SHBs llDa 



° This range includes the estimated thresholds of all the instru- 
ments that contributed to the detection of our SHB sample. 
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l2005|l . The fact that SHBs are observed in eUiptical 
galaxies cannot a priory rule out this possibility since it 
might be that neutron stars (and thus SGRs) can still be 
produced long after star formation ended in such galax- 
ies (Dar 20059- For example, it has been proposed that 
neutron stars can be the end produc t of white dwarf bi- 
nary mergers llSaio fc NomotolllQSl or of an accretion- 
induced collapse of a white dwarf ijNomoto fc Kondol 
I1991|) . However, observations indicate that SGRs and 
anomalous x-ray pulsars (AXPs) - two known popula- 
tions of magnetar candidates in our Galaxy - are asso - 
ciated with star forming regions l)Gaensler et a,l. ll2001|) . 
Thus our result that SHBs are dominated by progenitors 
with lifetime of several Gyr disfavor this possibility. 

3.1.2. Broken power-law luminosity Junctions 

Although a broken power-law luminosity function is 
not necessary in order to explain the data, this, or even 
more complicated functional forms, cannot be excluded. 
Given the prominence of broken power-law luminosity 
functions in previous studies we also consider a luminos- 
ity function of the form: 

m « { IZ\ llli^ (9) 

where ai < 1.5 and a2 > 2.5. We choose this range of 
indices so the luminosity function will significantly de- 
viate from the single power-law for which we already 
obtained a constraint on the lifetime. Another reason 
for exploring a broken power-law distribution is that it 
demonstrates the analysis that should be done once the 
sample of Swift short bursts with known redshift will be 
significantly larger and the Swift threshold will be bet- 
ter understood. Such luminosity functions couple the 
intrinsic functions to the detector response when Eq. ^ 
is integrated over L or z. Therefore, in this case, S{P) 
needs to be understood in order to compare the obser- 
vations with the model. Since we cannot construct S{P) 
for the whole sample we use in this section only the three 
bursts detected by Swift. However, as we now account 
for S{P) (and its uncertainties), we can use the two- 
dimensional distribution (Eq. Based on the peak 
flux of the three Swift SHBs with known redshifts we 
estimate the Swift/BAT threshold for localizations to be 
comparable to that of BATSE. We therefore approximate 
Swift's threshold as a step function at 1 ph/cm^/s with a 
large uncertainty due to the small number of Swift SHBs. 
We use Eq. 13 for the star formation history. 

In order to test the effect of broken power-law lumi- 
nosity functions on the derived limits on the lifetime we 
carry a maximal likelihood analysis for the narrow log- 
normal and the single power-law lifetime distributions. 
For each model we first fix ai = 1.5 and a2 = 2.5 and 
fit to the BATSE logN-logS distribution. A narrow 
lognormal lifetime distribution results with a most likely 
delay of = 6.5 Gyr and > 4[1] Gyr at 95%[99.98%]. 
A power-law lifetime distribution results with a most 
likely power-law index of 77 = 1 and with 77 > 0[— 1] 
at 95% [99.9%]. Repeating the analysis with other val- 
ues of ai and a2 (fitting to the BATSE logN-logS 
distribution) we obtain similar results. Given the un- 
certainty in the Swift threshold, we have also varied the 
assumed threshold value by a factor of 5 below our best 
estimate (retaining its step function form) , again finding 
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Fig. 4. — A color plot of the two-dimensional distribution 
dlog(L)dz ^^'^ '^^ "'^^ dimensional projections (obtained by inte- 
grating over z or log{L)) as predicted by a model with a lifetime 
distribution /(t) oc 1/t and an SFH of the form given by Eq. [7| 
The luminosity function is a broken power-law with «i = 1.5 and 
02 = 2.5 and L* = 2 ■ 10^^ erg/s, chosen to fit the BATSE LogN- 
LogS distribution. The threshold of Swift is taken as a step function 
at 1 ph/cm^/s. The Three Swift SHBs are marked with dots (on 
the two-dimensional plot) and bars (in the two projections). This 
model is rejected at 99.9% confidence level. 
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T.=3 Gyr ; a=0.3 

L.=10^^ erg/s 
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a2=2.5 
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Fig. 5. — Same as Fig. |^for a model with a lognormal lifetime 
distribution (r« = 3 Gyr and a = 0.3) and L» = lO^'^ erg/s (the 
rest of the parameters are the same as in Fig. |3J. This model is 
rejected at 99% confidence level. 

similar results. Note that although here we can use a 
significantly smaller sample (3 bursts) than the one we 
used in H3.1.1I (8 bursts), we obtain similar significance. 
The reason is, at least partially, the usage of the two- 
dimensional distribution which contains more informa- 
tion than the redshift distribution alone. This sample is 
however much more vulnerable to unaccounted errors^. 

Figures 0] and [SI depicts the predicted two-dimensional 
distributions of a power-law lifetime distribution with 

— ~1 and a narrow lognormal lifetime distribution with 
— 3 Gyr, respectively. These figures include also the 
observed distribution as well as the two one-dimensional 

^ Repeating the analysis when the redshift of SHB 0508013 is 
taken as a lower limit > 0.72 (instead of = 0.72) does not alter 
the main results, that a long lifetime is required. The most likely 
values of t, and rj are similar while r* < 1 Gyr is rejected at 99.7% 
and r; < — 1 is rejected at 99.4% 
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t:=6 Gyr ; a=0.3 

L.=4- 10^° erg/s 

a^=1.5 

a2=2.5 
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Fig. 6. — Same as Fig. |3 for a model with a lognormal lifetime 
distribution (t« = 6 Gyr and a = 0.3) and L» = 4 ■ 10^*^ erg/s 
(the rest of the parameters are the same as in Fig. EJ. This model 
provides a good description of the data. 

projections of the models and the observations. The fig- 
ures show clearly that even though there are only three 
observed SHBs in this sample the two models do not de- 
scribe the data well. Indeed both models are rejected at 
a confidence level that is larger than 99%. 

This result can be understood. Given the intrinsic 
rate evolution Rshb{z) predicted by lifetime distribu- 
tions with short delay (< 3 Gyr) or small power-law in- 
dex (77 < —0.5), the luminosity function (Eq. (S)) implies 
that the typical observed burst luminosity is L*, and 
the typical observed redshift is the redshift at which 
is detected at the threshold level^". Regardless of the 
exact value of ai and a2 we find a typical lumi nosity 
> 10^^. Similar values of were obtained bv I Andol 
|2004) for various luminosity functions. The typical red- 
shift is z « 0.7, and together with these two values 
are too high for the Swift sample. Longer lifetimes result 
in Rshb{z) that increases at low redshifts, leading to a 
lower value of as well as to lower typical redshifts. For 
example, a delay of 6 Gyr results in = 4 x 10^° erg/s 
and predicts a distribution that is fully consistent with 
the observations (see fig. EJ. 

Finally we discuss the robustness of our limit on the 
progenitor lifetime in the case of a luminosity function of 
an arbitrary form. The reason that short lifetime models 
fail is the deficit of low luminosity SHBs at low redshifts. 
In order to increase this fraction the luminosity function 
must increase more steeply than at low luminosities 
(< 10^0 erg/s). In this case however the luminosity func- 
tion must become fiatter than L^^ at higher luminosities 
in order to fit the BATSE LogN-logS (i.e the luminosity 
function must have an 'ankle'). Although we cannot ex- 
clude a luminosity function with an 'ankle', and thus the 
possibility of a shorter typical lifetime, such a luminosity 
function is not expected unless SHBs are composed of 
two separate populations. 

An alternative method to estimate the progenitor life- 
tim e is based on t he spectral types of the host galax- 
ies (|Gal-Yam e,t al. ..2005') . In the extended sample four 
out of the six putative host galaxies are of early type 

This statement is correct as long as L* is detected up to a 
redshift < 1, which is the case for the values of L, that we find 
and the Swift threshold. 



If only Swift and HETE2 bursts are considered then 
three out of four are of early-type. The two late-type 
galaxies are moderately star- forming In one a signif- 
icant fraction of th e stars is ^ 1 Gyr old (SHB 050709; 
iCovino et al.ll2f)r)6h . In the other most of the stars are 
much older than 1 Gyr (S HB 000607: IGal-Yam et al. I 
12005ft . Using these results IcTl-Yam et al. I l|2005j) find 
that SHBs are more likely to reside in early type galax- 
ies than type la supernovae (at 2a c.l.), implying that the 
typical lifetime of SHB progenitors is most likely longer 
than the one of type la SNe (> 1 Gyr). This method 
does not depend on the luminosity function or the instru- 
mental thresholds and thus it constitutes an independent 
corroboration of our results. 

3.2. The Local Rate of SHBs 

The most robust lower limit on the local SHB rate is 
obtained directly from the BATSE observed rate (a full 
sky rate of ~ 170 yr~^). Considering a single power-law 
luminosity function the observed redshift distribution of 
the extended sample (8 bursts) should be similar to the 
BATSE redshift distribution of BATSE bursts. There- 
fore, the fact that 5 out of the 8 bursts are within a 
distance of 1 Gpc suggests that the rate of SHBs ob- 
served by BATSE is « 25 Gpc~^ yr~^ and puts a lower 
limit of ~ 12[6] Gpc^^ yr~^ at 2[3]a confidence level. 
Even with no assumptions about the luminosity func- 
tion, the sample of SHBs detected by Swift (3 bursts), 
which has a comparable threshold to that of BATSE, 
gives a similar result. Two of these bursts are within 
1 Gpc suggesting the same rate of ~ 25 Gpc~^ yr~^ 
and setting a lower limit of « 6[1] Gpc"'^ yr~^ at 2[3]cr 
confidence level. The rate of SHBs observed by BATSE 
is a strict lower limit since it does not include unde- 
tected SHBs. Bursts can avoid detection either by point- 
ing away from us, if the prompt gamma-ray emission is 
beamed, or by being too dim. This robust lower limit 
is higher than previous estimates (fSchmidtll2001t lAndol 
2004; Guctta & Piran 2005) by a factor of 10-100. This 
observed local rate is also higher by the sa me factor than 
the observed local rate of long-so ft GRBs ljSchmidtlll999t 
IGuetta,. Piran Wa,xma,nll200l. 

The progenitors of SHBs are almost certainly a product 
of at least one core-collapse supernova (e.g., a neutron 
star) and the SHB itself is most likely a catastrophic 
event, therefore the rate of SHBs is limited by the rate of 
core-collapse SNe. Since we find that the typical lifetime 
of SHB progenitors is several Gyr the local rate of SHBs 
is limi ted by the rate of co re collapse SNe at a redshift 
« 0.7 ijPahlen et al. Il2004j) . Together with the observed 

Note that the host of SHB 050813 is of early-type re gardless of 
whether it is associated with the cluster at z=0.72 iGladders et g^l 
12005; Bcreer 2005; Prochas ka et al. 120051) or with a cluster at 2 Ri 
1.8 (Bcrgcr 2006). Moreover, if the burst is truly associated with a 
cluster at 2 Si 1.8 than the lifetime of its progenitor is 3 Gyr. It 
cannot be older since this is the age of the universe at this redshift 
and it is not much younger since it is associated with one of the 
most developed, and therefore one of the oldest, structures at this 
redshift. 

Note that as we discuss in JSthe host of SHB 051221 cannot 
be included in this sample without accounting for selection effects 
that prevent the detection of the hosts of other SHBs that were 
localized after we defined our sample. This is especially important 
when the host type is considered since we expect afterglows to be 
bright in a gas rich environments, making the identification of a 
late-type hosts more probable. 
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lower limit we find that the local rate of SHBs, TZshb = 
Rshb{z = 0), is in the range: 

10 < 7^sHs < 5 X 10^ Gpc"^ yr~^ (10) 

Constraining the local rate within this range requires 
an estimate of the beaming factor (i.e., the fraction of 
the 47r solid angle into which the prompt gamma-rays are 
emitted), f^^^ and of the luminosity function. The value 
of the beaming factor is currently unknown, however, 
the afterglows of two SHBs (050709 and 050724) have 
shown a steepening that can be interpreted as a hint of 
a jet (Fox et al. 2005). This interpretation, although 
somewhat speculative at this point, indicates a beaming 
factor of ^ 50. The correction for undetected dim bursts 
depends on the luminosity function and most strongly on 
its lower cutoff, Lmm- Taking the power-law luminosity 
function that we have found above {(j>{L) oc i~^) the 
local rate is weakly sensitive to the exact evolution of 
RsHB with redshift (as long as the fraction of nearby 
bursts is compatible with the current sample ^■^): 

Taking the beaming factor suggested by iFox et aL\ 
(j2005) (30-50) and the value of Lmi„ sug gested by our 
results in conjunction with the results of iTanvir et al. I 
p)05li) (lO^^erg/s), we obtain a rate of 10^ Gpc'^ yr'^. 
This high rate is highly uncertain at this point but we 
consider it to be the most favored by the synthesis of cur- 
rent observations. This rate is only slightly smaller than 
the upper limit we obtained (Eq. 1101) . implying that if 
this estimate is correct a significant fraction of all the 
end products of core-collapse SNe are producing SHBs. 
Alternatively, perhaps a source can produce more than a 
single SHB (i.e. SHBs are not catastrophic events). The 
upper limit that we find does not pose at the moment a 
stringent limit on the beaming or on Lmin- However, if 
in the future bursts with a luminosity of 10^^ erg/s will 
be observed, it will constrain the beaming to be smaller 
than ~ 100. 

4. SHBS AND COMPACT BINARY MERGERS 

4.1. The predicted lifetime and merger rate of DNS and 
BH-NS binaries 

The merger rates of DNS and BH-NS systems are 
currently estimated in two ways, based on observed 
systems in our galaxy and using theoretical population 
synthe sis. Esti mates based on observed DNS b i naries 
Ce.g.. iPhinnevI IT9 91: Naravai TTiran fc Sheiiiil 119911: 
Curran & Lorimcr 1995: van den Hcuvel & Lorimer 
19961 lAr zou mania n. Cordes & Wasscrman 199E; 
Kalogera, et nH [20011 I200I Ide Freitas Pa,checo et aU 
2005j) provide a lower limit on the local rate of DNS 
merger events because of possible unaccounted for selec- 
tion affects (e.g., binary formation that does not recycle 
neither of the neutron stars). The absence of observed 
BH-NS systems together with poor understanding of 
selection affects involved in their discovery renders this 

The rate varies by a: 30% when calculated using the most 
likely values of the three models that we consider. Adopting values 
values of the models that are rejected at 2[3]cr can reduce the rate 
by a factor of 3 [6]. 



method u napplicable in this case. Population syntheses 
(e,g. Linunov et al 111 99,'^ IPortegies Zwart fc Yungelson 
19981 iBethe fc Brownl 119981 iBloom. Sigurdsson fc Pols 
iFrve r^'^oos leTir^artmam 

^9^ "Bclczvriski fc Kiifogeral 

Belczvnski. Ka logcra fc BuliB |2QM 

Belczvnski. Buli k fc Kalogeral l2002t iPerna fc B elczvnskil 
do not suffer from observational selection effects 
and can address both DNS and BH-NS binaries. How- 
ever, the uncertainties involved are substantial, with 
rate estimates spanning over three orders of magnitude. 

The estimates based on observed DNS systems were 
drastically revised with the discovery of the relativis - 
tic binary pulsar PSR J0737-3039 (Bur gav et aZ. 120031). 
Based on three observed DNS binaries^^ iKalogera et al. I 
(|M04) find a range of Galactic merger rates from 1.7 x 
10-5 to 2.9 X lO-'^yr-i at 95% confidence. This es- 
timate is larger by a factor of 6-7 than the estimates 
that they obtained excluding PSR J0737-3039. Extrap- 
olation of this rate to the cosmological neighborhood 
(jKalogera et al. , 20M) yields LIGO-I and LIGO-II de- 
tection rates of 7 x 10"^ - 0.1 yr'^ and 40 - 650 yr"^ , 
respectively. This rate corresponds to a local rate density 
of 200- 3000 Gpc-3 yr-i. 

Regardless of the exact values of the above esti- 
mated rate, the merger rate for observed DNS binaries 
is dominated by systems with short lifetimes. This is 
demonstrated, for example, by the fact that the detec- 
tion of PSR J0737— 3039 (with a merger timescale of 
« 100 Myr) significantly increased the estimated merger 
rate (Kalogcra et al. 2004). Note that this drastic re- 
vision in the rate is partially due to the short lifetime 
of this double pulsar. Numerous new binaries with long 
lifetimes need to be discovered if old binaries are to dom- 
inate the merger rate in our galaxy. 

The rates found by population syntheses for DNS sys- 
tems are in general consistent with the rate deduced from 
the observations. Analysis of BH-NS binaries shows that 
their merger rate may be comparable, or even larger by 
an order of magnitude l)Bethe fc Brownl ll99^. than the 
rate of DNS mergers. Population synthesis provides a 
rather robust upper limit on the rate of DNS and BH- 
NS mergers by considering the fraction of binaries that 
survive two core-c ollaps e SNe. Several theoretical works 
(e.g^lPfahl et al. Il200l iLipunov. Postnov fc Prokhorovl 
Il997^ show that this fraction is likely to be ~ 0.001 and 
can be as high as 0.02. Tak ing this fraction from the local 
rate of core-collapse SNe ijCappellaro. Evans fc Turattd 
II999I) implies an upper limit of 1000 Gpc"'^ yr~^ on the 
merger rate of short-lived binaries (t < 1 Gyr). An upper 
limit of 10^ Gpc~'^ yr"^ on the merger rate of long-lived 
binaries (t > 4 Gyr) is obtained by considering the rate 
of core-collapse SNe at redshift 0.7 (jDahlen et al. 1.2004^) . 

4.2. Comparison with SHBs 

DNS or BH-NS mergers are currently the most 
popular models for SHBs ( e.^.. lEichler et al. I 119891 
iNaravan. Paczvnski fc PiranI 119921. The proposed 
association is based on the tight limits that are imposed 
on the progenitor by the large energy release 10^° 

IKalogera et al.^ 120041) consider only the observed DNS sys- 
tems that will merge within a Hubble time. They further exclude 
PSR B2127-I-11C because of its association with a globular cluster. 



erg) and the short timescales 1 nis) observed in 
these bursts. Not many sources can produce such an 
electromagnetic d isplay, while an exte n sive numeri- 
cal e ffort Ce.g. iRuffert fc Jankal 119991: iJanka et aL") 
19991 iRoss wog. Ramirez-Ruiz fc Davics" '2003^ 
Lee. R,amirez-R,uiz fc Granot 20Q5: Oechslin fc Janka 
2005ft has shown that the mergers of DNS and BH-NS 
systems can. 

The range of SHB rates that we find (Eq. Ill|) overlaps 
the predicted rate of DNS mergers based on observations 
(~ 200 — 3000 Gpc~^ yr~^). However, the range of pro- 
genitor lifetimes we find for SHBs is inconsistent with 
that of the observed DNS binaries. The reason is that 
the merger rate of the latter is dominated by short-lived 
binaries (r 100 Myr) while SHB progenitors are much 
older (r > 4Gyr). We do not see any obvious selection 
effect that may operate and can be applied to the cur- 
rent DNS models and explain this discrepancy. On the 
contrary, while very tight DNS binaries (with very short 
lifetimes) can avoid detection in pulsar surveys because 
of their large orbital accelerations, long-lived systems are 
detected with relative ease as long as one of the neutron 
stars is an active pulsar. 

Admittedly, the probability that both members of a 
DNS binary are dead pulsars (i.e., undetectable), in- 
creases with its age. However, our limits on the life- 
time distribution of SHBs constrain it to be shallower 
than r"°'^ (at 95% confidence level). In other words, 
the number of newly born SHB progenitors with a life- 
time between 1 — 10 Gyr is larger by a factor > 10 than 
those with hfetime 10-100 Myr. Therefore, if DNS bina- 
ries are the progenitors of SHBs, for every recently born 
DNS with a lifetime 10 — 100 Myr we expect to observe 
at least several recently born systems with a lifetime of 
1—10 Gyr. All of these coeval systems are presumably 
equally detectable via pulsar surveys. In reality there 
is one observed DNS with a lifetime of 10 — 100 Myr 
and only one with a lifetime of 1 — 10 Gyr (B1534-I-12, 
which is most likely older than 100 Myr). We therefore 
conclude that with the caveat of small number statistics, 
the observed population of DNS systems is unlikely to 
be the dominant component among SHB progenitors. 

It is possible that known DNS pulsars are not rep- 
resentative of the larger population of relativistic DNS 
systems. This might be the case if large numbers of 
relativistic DNS systems are formed without recycled 
pulsars, or alternatively if the Milky-way DNS popula- 
tion do not represent the cosmic one. We thus conclude 
that SHB observations disfavor a DNS merger origin to 
the extent that they indicate that an unexpected hidden 
population of old long-lived systems should be invoked. 
While we cannot exclude the existence of such a pop- 
ulation at the moment, it remains to be seen whether 
it can be reproduced by population synthesis models. 
Turning the argument around, if DNS binaries are the 
progenitors of SHBs then there exists a large popula- 
tion of hitherto-undetected old DNS systems. More- 
over, this scenario suggests that the rate of SHBs is 
TZ ~ 10"^ — lO'' Gpc^'^ yr^^ in order to dominate over 
the the merger rate o f young DNS systems found by 
iKalogera et al. I l)2flfl4D . while conforming with the up- 
per limit imposed by the rate of core-collapse SNe at 
2 « 0.7 and the binary survival probability. It also ex- 
cludes the existence of the recently-hypothesized dom- 
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inant population of DNS binaries with very short life- 
times (|Belczvnski fc Kalogerall2001t iPerna fc Belczvnskil 
I2OOI . 

The possibility that BH-NS systems are the dominant 
source of SHBs cannot be constrained at this time. There 
are no observational data on the properties of this popu- 
lation while the theoretical models are highly uncertain 
and their predictions vary significantly from one model 
to the next. If a link between SHBs and BH-NS merg- 
ers will be established in the future, the observed SHB 
properties can be used to constrain models of NS-BH 
binaries. 

5. PROSPECTS FOR GRAVITATIONAL WAVE DETECTION 

Progenitors of SHBs may also be strong sources of 
gravitational waves (GWs). Therefore, the local SHB 
rate might has direct implications for the detection rate 
of GW telescopes in general and to the upcoming S5 
run of initial LIGO (LIGO-I) in particular. The most 
promising progenitors for GW detection are NS-BH or 
DNS mergers. As we discuss above our results disfavor, 
but cannot exclude, DNS binaries as the source of SHBs, 
while the possibility that NS-BH binaries are the progen- 
itors on SHBs cannot be constrained by the data. More- 
over, energy and time scale considerations favor compact 
binary mergers as the progenitors of SHBs. Therefore it 
is worthwhile to discuss the predictions for the LIGO 
and VIRGO GW observatories if SHBs originates during 
such compact mergers. 

These sources are expected to be detected by LIGO-I 
up to a distance of - 43 [20] Mpc ijCutler fc Thornd2002t 
hereafter the numbers without brackets are for NS-BH 
mergers assuming a BH mass of lOM© while the numbers 
within the brackets are for DNS mergers^^). VIRGO de- 
sign sensitivity is similar to that of LIGO-I, and there- 
fore the VIRGO detection rate is expected to be similar 
once its design sensitivity is achieved. If SHBs are trig- 
gered by compact mergers, the prospects for GW detec- 
tion from SHBs are quite promising. A reasonable, but 
speculative, SHB rate of 10** Gpc"^ yr"* (see 
predicts a detection rate of ~ 3 [0.3] mergers per year. 
This speculative rate is higher by a factor of ~ 10 than 
previous upper limits on merger rates derived from the 
local rate of core-collapse SNe. The reason is that the 
long lifetime of SHB progenitors relates them to SNe at 
redshift ^ 1, which was higher by a factor of ^ 10 than 
it is today. 

The possible temporal coincidence between an SHB 
and a GW signal, which may make the difference be- 
tween detection and non-detection in the coming S5 run 
of LIGO-I. A coincident detection increases the range 
of LIGO-I and VIRGO by a factor of w 1.5, while if 
SHBs are preferentially beamed perpendicular to the bi- 
nary orbital plane the range is increased by an addi- 
tional factor of » 1.5 to a distance of « 100 [50] Mpc 
l|Kochanek fc PiranllT99^ . Hence LIGO-I has a non- 
negligible probability to detect GWs simultaneously with 

15 Recently iMiiieJ 120(1^ and lR,osswo3 1200.'^^ argued that BH- 
NS mergers with a mass ratio > 10 cannot produce the accretion 
disk that is believed to be necessary for the production of SHBs. 
This result suggests that SHBs can be produced only by mergers of 
BH and NS with similar masses. In this case the W detectability 
of BH-NS mergers is similar to that of DNS mergers, 
http:/ /www. virgo.infn.it/senscurve/ 
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a Swift detection of an SHB. Our best fit model to the 
current observations {(p{L) cx L^^, Lmin < 10^'' erg/s 
and any lifetime distribution that provides a large frac- 
tion of SHBs at low redshift) predicts that about 3% of 
the SHBs observed by BATSE are within 100 Mpc while 
about 1% are within 50 Mpc. Given that Swift detects 
^ 10 SHBs a year and that its threshold is comparable 
to that of BATSE, the expected rate of simultaneous de- 
tections of an SHB and a GW signal is ~ 0.3[0.1] yr~^. 
This rate depends on the beaming only through the as- 
sumption that SHBs are beamed perpendicular to the 
binary orbital plane. It could be further increased by a 
search for simultaneous GW signals and SHBs that are 
detected by Swift or the Interplanetary Network (IPN) 
but are not localized (such bursts are currently not an- 
nounced and are not included in the Swift detection rate 
that we quote above). 

Simultaneous detection of the inspiral GW signal from 
a compact merger and an SHB will provide a conclusive 
evidence that SHBs originate from compact mergers and 
would improve our understanding of both merger physics 
and SHBs significantly. GWs can provide a unique view 
of the formation, and possibly the operation, of the in- 
ner engine powering burst, difficult to observe via any 
other observational method. If GWs from a merger are 
detected without a detected SHB, an association may 
still be se cured by the detectio n of an on-axis orphan 
afterglow l|Nakar fc Piranll2003|) . The probability to de- 
tect an off-axis radio orphan afterglow as suggested by 
iLevinson et al. I l|20 02') is low unless the beaming is not 
significant. The reason is that the isotropic energy of an 
SHB that would be detected by LIGO-I or VIRGO is ex- 
pected to be < 10**^ erg, and if the beaming is significant, 
its total energy would be much lower, and the afterglow 
is expected to be too dim for detection by the time that 
the jet decelerates and spherical emission dominates. 

Advanced LIGO (LIGO-II) would be able to detect 
NS-BH [DNS] mer gers up to a distance of 650 [300] Mpc 
ijCutler fc Thornd 12002). Thus, if SHBs originate in 
such mergers then a robust lower limit for the detec- 
tion of GWs by LIGO-II is 20[2] yr^^ and the probable 
rate is larger by orders of magnitude. This rate is for 
detections of GW signals that are not associated with 
the prompt gamma-ray emission from an SHB. As dis- 
cussed above, simultaneous detections will increase the 
LIGO-II range by a facto r of 1.5-2.5 to « 1.3[0.6] Gpc 
l|Kochanek fc Piranll99.? ). So far, T/E^TE^-g observed one 
burst at a distance < 700 Mpc while Swift detected at 
least additional 2 SHBs at distance < 1 Gpc. The GBM 
detector on GLAST^'' is expected to have a threshold 
that is similar to BATSE and more than half sky field of 
view, and thus it is expected to detect at least 5 SHBs 
within a distance of 500 Mpc every year. The benefit of 
the simultaneous operation of LIGO-II and an efficient 
GRB detector goes beyond the high likelihood to observe 
simultaneous SHBs and mergers if they are associated - 
it will also be able to disprove the association at high 
significance if no simultaneous detection is observed. 

SHBs may emit GWs also if the progenitor is not a 
compact binary merger. Any progenitor that involves 
the collapse of a rotating compact object to a black hole 
(e.g., the collapse of a rotating neutron star triggered 

'-^ http:/ /f64. nsstc.nasa.gov/gbm/ 



by accretion; iMacFadven. Ramirez-Ruiz fc Zha ndl2005{) 
will p roduce gravitational waves fsee. e.g.. lStark fc PiranI 
Il985|) . The amplitude of these waves is highly uncertain. 
Moreover, the absence of an accurate signal template will 
reduce its detectability. Such GW signals would most 
likely not be detected even by LIGO-II at distances m uch 
greater than 10 Mpc ijKokkotas fc Stergioulagll2'005l and 
references therein). If SHBs result from the collapse of a 
compact object that is not triggered by a merger, LIGO- 
II might detect several SHBs if they are significantly 
beamed and/or Lmin <C lO''^ erg/s. 

6. CONCLUSIONS 

We have extended previous analysis schemes (jPiranI 
119921 lAndo 2004: Guetta fc Piran 2005), combining the 
observed peak flux distribution and the two-dimensional 
observed redshift-luminosity distribution of SHBs, to 
constrain their progenitor lifetime and their local rate. 
We apply this method to the extended sample of SHBs 
with redshift constraints presented in iGal-Yam et al. I 
12005). Our conclusions are: 

• The progenitors of SHBs are dominated by an old 
population. If there is a typical lifetime than it 
is > 4[1] Gyr at 95%[99.9%] confidence level. If 
the lifetime distribution is a power- law with index 77 
then the number of progenitors per logarithmic life- 
time interval increases significantly with lifetime - 
rj > -0.5[-l] at 95%[99.5%] confidence level. Sim- 
ilar results are obtained when only the Swift SHB 
sample is considered. These results were obtained 
assuming that the luminosity function and the life- 
time distribution arc unimodal. 

• We derive a lower limit on the local rate of SHBs, 
T^SHB.obs ^ 10 Gpc^^ yr^^- This rate is compa- 
rable to the BATSE observed (all-sky) rate and is 
a robust lower limit. It is obtained when we use 
the simplest luminosity function that fits the data, 
(j)L oc L~^, or independently, if we take the Swift 
threshold to be comparable to that of BATSE, as 
implied by the SHBs observed by Swift so far, and 
consider only the Swift sample. 

• Assuming that SHB progenitors are the end prod- 
ucts of core-collapse SNe, and that they are catas- 
trophic, we derive also a rate upper limit - the rate 
of core-collapse SNe at z w 0.7: 10 < TZshb < 
5 X 10^ Gpc-3yr-i. 

• Indications that SHBs are beamed l|Fox et al. I 
12005^) and that the luminosity funct ion lower cut- 
off is Lmin < 10''^ erg/s f^^O and iTanvir eit al I 
()2005D ) suggests that the local rate is much higher 
than our lower limit. Considering the current 
best, but highly uncertain, estimates of the beam- 
ing and Lmin we get a local rate of TZshb ~ 
10^ Gpc-3yr-i. 

We compare our constraints on the progenitors of SHBs 
to observational constraints on DNS mergers. We find 
that observed DNS binaries in our galaxy are unlikely 
to be a representative sample of SHB progenitors. The 
reason is that the merger rate of the observed DNS bina- 
ries is dominated by short-lived systems (r < 100 Myr) 
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while SHB progenitors are much older (r > 4Gyr). We 
are not aware of any obvious selection effect that would 
prevent the detection of recently formed yet long-lived 
DNS binaries and therefore we conclude that our results 
disfavor DNS systems as the progenitors of SHBs. How- 
ever, given the limited understanding of the formation of 
DNS systems and the possibility that some of the forma- 
tion channels result in binaries that cannot be detected, 
we cannot exclude the possibility that DNS systems end 
their life as SHBs. If they do, then we expect a large 
number of undetected old DNS systems in our galaxy. 
We consider also the theoretical constraints on BH-NS 
mergers and find that there are no robust constraints 
that prevent these from being the progenitors of SHBs 
(based on the lifetime and the rate) , due to the uncertain 
properties of the BH-NS population. 

If SHBs originate from either BH-NS or DNS mergers, 
then the prospects for detecting GWs by first genera- 
tion detectors and in the future planned advance LIGO 
(LIGO-II) are promising: 

• Our robust lower limit on the local rate of SHBs 
implies that LIGO-II will detect at least 20 merg- 
ing BH and NS pairs with a mass ratio of a few, 
per year, if they are the progenitors of SHBs. A 
minimal rate of 2 DNS (or BH-NS with similar 
masses) mergers per year is predicted if these are 
the sources of SHBs. If an efficient GRB detector 
(e.g. GLAST) will be operational contemporane- 
ously with LIGO-II, numerous simultaneous detec- 
tions are expected. The absence of such detections 
will disprove models of SHBs are as compact binary 
mergers. 

• The highest possible rate of SHBs that is consistent 
with the upper limit on the merger rate of compact 



binaries (10"* Gyr yr ^) predicts a LIGO-I and 
VIRGO detection rate of - 3 yv~^ if SHBs are 
mergers of BH-NS with a mass ratio of a few, and 
~ 0.3 yr~^ if they are DNS mergers or mergers 
of BH-NS with similar masses. This is the rate 
for GW detection only, with no associated SHB 
detection by Swift. 

• The predicted probability for a simultaneous detec- 
tion of GWs by LIGO-I or VIRGO and an SHB by 
Swift in one year of contemporary operation might 
not be negligible: '-^ 30% for a BH-NS merger and 
~ 10% for a DNS merger. This probability is de- 
rived assuming (f>L oc with a lower luminosity 
cutoff L„^in = 10'*'^ erg/s. 

• The probability for a simultaneous detection can 
be increased if a search for GW signals will be car- 
ried out at times that SHBs are detected, but not 
localized, by Swift or by any of the IPN spacecraft. 
If SHBs are not strongly beamed, this may signifi- 
cantly increase the probability of LIGO-I to detect 
a single event. 
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